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Feldbacher, Held, and Assaad Reply: Re- 
cently, we introduced a projective quantum Monte Carlo 
(PQMC) algorithm for simulating the Anderson impurity 
model (AIM) Q . The preceding critique [2| based on the 
orthogonality catastrophe (OC) Q is invalid: (i) There 
is no OC in 0], and it is generally not "unpractical" to 
avoid it. (ii) The OC does not affect our results. 

Concerning (i): The OC theorem [Sj states: If a Hamil- 
tonian Ht is perturbed by a local disorder potential 
and/or interaction H' to a Hamiltonian H = Ht + H' , 
the overlap of the ground states of H and Ht is 

KV'gIV't)! ~ 7V~",7Vbeing the numer of bath sites. (1) 

What Katsnelson 2| overlooked is that a = in pl . The 
a of Eq. Q is given by the difference between the scatter- 
ing phase for Ht and H (tpT and ipa, respectively), i.e., 
a ~ {ipT — Vg)^|J|- According to the Friedel sum rule, 
the scattering phase ip is related to the average number 
of electrons on the impurity site: f — tt/2 n^. Since n^ 
is the same for H and H' , we have (pT — ^g, a-nd hence 
a = 0. There is no OC. 

Note that our calculations in [l|, |5j are for half-filled 
bands, so that we automatically have a = 0. But also off 
half- filling one is free to choose a Ht with the same n^j. 

This proof based on Friedel's sum rule has been incor- 
porated in |2|, where it is now argued that using a Ht 
without OC is "unpractical" |2| because one would need 
an (analytically) "exact answer for n^" 2\. We reject 
this critique. For a numerical algorithm, it is perfectly 
legitimate to calculate rid numerically. We have done so 
in practice j6j. 

Concerning (ii): The OC is irrelevant for our calcula- 
tions since for long enough projection time 9 we obtain 
the same Green function G{t) with and without OC, see 
Fig. 1. The calculations of Fig. 1 have been done for an 
AIM with on-site hybridization F = 1, a constant bath 
density of states ranging from -1 to 1 (with A^ — > oo bath 
sites), imaginary time discretization Ar = 0.1, impurity 
level ed = 0.5, and Coulomb interaction U — (for maxi- 
mizing precision). The trial AIM with OC was half-filled, 
i.e., e*/^! = 0, n*™i = 1 instead of e^ = 0.5, Ud ~ 0.7. For 
long enough projection 0, the Green function with OC 
converges to the exact result. At 6 = 100 (250), the av- 
erage difference to the exact G{t) is already as small as 
7-10-'' (3-10-''). Also, rid = 2[1-G(t = 0)] is obtained 
extremely precisely even with OC. 

Technically, one can understand the irrelevance of the 
OC by considering that we have integrated out the N 
bath sites of Eq. Q [ll l3| . The algorithm is exactly the 
same with OC {N — > oo) and without OC (finite N). For 
large enough but finite N, also the non-interacting input 
Green function can be the same as for N -^ oo, within 
any given accuracy. Hence, the results are the same with 
and without OC. 

Physically, one can understand the irrelevance of the 
OC by considering that the OC theorem 3^ holds in the 
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FIG. 1: (Color online) PQMC Green function with OC and 
without OC (exact). For long enough projection 0, the Green 
function with OC converges to the exact result. Note the 
logarithmic scale of the y axis. 


metallic, (quasi-)particle regime and only for N — > oo. 
Then, however, the ground state and excited states with 
a finite number of particle- hole excitations, have the same 
physical properties, since the finite number of particle- 
hole excitations becomes irrelevant for N —f oo. Ground 
state and such low-lying excited states yield the same 
equilibrium (ground state) Green function Q. This ex- 
plains, why we obtain the correct results, e.g., in Fig. 1, 
even with OC. The OC is irrelevant for our calculations. 
Katsnelson's objection J3| is not valid. 
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